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ABSTRACT 


This document preijents the mathematical formulation of CNTRLF, the maneu- 
ver control program for the Applications Technology Satellite-F (ATS-F). The 
purpose of the document is to specify the mathematical models tkat are included 
in che design of CNTRLF, It serves as final documentation, superseding the 
preliminary Mathematical and Programming Task Specifications for the ATS-F 
Orbital Maneuver Control Program (CNTRLF). 
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SECTION 1 - INTRODUCTION 


This document specifies the mathematical models that are included in the de- 
sign of CNTRLE — the maneuver control program for the Applications Technol- 
ogy Satellite-6 (ATS-6). ^ Also, the ATS- 6 maneuver control hardware and 
maneuver requirements are summarized to identify the control system char- 
acteristics which were modeled. 

To satisfy the specific requirements of different audiences, CNTRLF is 
documented in several volumes. The Mathematloal and Programming Task 
Specification for the ATS-F Orbital Maneuver Control Program (CNTRLF) 
(Reference 1) describes the preliminary functional requirements and design 
of CNTRLF. It provides the analyst and programmer with the background 
material needed to understand the evolution of CNTRLF. The Subsystem 
Description for the Applications Technology Satellite- F (ATS- F) Orbital Ma- 
neuver Control Program (CNTRLF) (Reference 2) contains the details of the 
CNTRLF design and is directed toward programming and maintenance person- 
nel. The User’s Manual for the ATS-F Orbital Maneuver Control Program 
( CNTRLF) (Reference 3) is directed toward a general audience of analysts, 
programmers, and data technicians. In addition to a brief description of the 
program, the principal contents of the User’s Manual are descriptions of the 
specific requirements for data card input to the program. 

This document is final documentation and supersedes Reference 1. It serves 
as an introduction to CNTRLF for programmers who must maintain and en- 
hance the CNTRLF program. It is most helpful also for the analyst who must 
be familiar with the program at the algorithm level. 


^ATS-F was renamed ATS-6 after its launch in May 1974. For consistency, 
ATS-6 will be used throughout the remainder of the document. 
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1.1 HISTORICAL BACKGROtnSiD 


Originally, design of the orbital maneuver control program for ATS-6 was to 
be based on the existing program for ATS series spacecraft (CONTRL)* Ref- 
erence 1 discusses this design based on modifications to the CONTRL program. 
The principal modifications to CONTRL included; the calculation of the thrust 
vector direction; the calculation of the change in magnitude of the velocity 
vector (AV) for drift rate and inclination corrections; the calculation of the 
propulsion system operating state during bums; the prediction of orbital 
changes resulting from ion engine operation; and the elimination of computa- 
tions required for attitude maneuvers. Early in the implementation of this 
design it became evident that the ATS-6 control problem was much more unique 
than originally envisioned. The scope of the modifications required to accom- 
modate the difference between the stabilijjation and propulsion subsystem of 
ATS-6 and the operational ATS-1 and -3 warranted changes to the original 
design. 

In a parallel effort, the structure of the Synchronous Meteorological Satellite 
(SMS) maneuver control program (SMSMAN) (Reference 4) was being developed, 
SMSMAN, a generalized maneuver control program for spin-stabilized satel- 
lites contained many of the features required in the design of CNTRLF. There- 
fore, rather than modifying CONTRL to support ATS-6, the SMSMAN program 
was chosen as the base program. 

The CNTRLF program is being used operationally to compute orbital adjust 
maneuvers for ATS-6. The original design has been updated to include station- 
keeping logic, node change maneuver logic, an improved method for calculating 
the drift rate, a hydrazine engine propulsion model, a finite bum integration 
model, and an improved internal orbit generator. 
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1.2 ATS-6 SPACECRAFT DESIGN FOR PERFORMING ORBITAL ADJUST 
MANEUVERS 

Orbital adjust maneuvers as well as certain attitude control functionB for 
AT3-6 are performed with the Spacecraft Propulsion System (SPS). Two 
experimental ion engines are available for supplementary inclination and node 
corrections. The SPS uses a hydrazine (N H .) thrusting system, consisting 

A 4 

of 2 fuel tanks pressurized by gaseous nitrogen (N ) and 16 heated thrusters. 
The SPS is completely redundant, and fuel from either tank may be directed 
to any nozzle by switching the latching valves incorporated in the fuel lines. 

The location and thrusting direction, with respect to the spacecraft's X , Y , 
and Z axes, of each of the SPS thrusters is known. Figure 1-1 shows the con- 
figuration of the thrusters on the orbital control jets (OCJ) bar. Thrusters 7, 

8, 15, and 16 are dedicated solely to east-west orbital maneuvering. Thrust- 
ers 13 and 5 and thrusters 14 and 6 are used in pairs to control the north-south 
orbital maneuvering. 

The two ion engines are mounted with their thrust axes nominally at ±45 de- 
grees from the spacecraft's +Z-axis in the Y-Z plane as illustrated in Fig- 
ure 1-2. In general, the thrusters are positioned so that their thrust acts 
through the center of gravity. However, some torque will be placed on the 
spacecraft due to changing foel levels and initial misalignments. 

1.2.1 Orbital Maneuvers With the SPS 

The SPS provides thrust to change the orbital velocity of the spacecraft for 
east-west and north- south orbit control. The SPS also provides thrust to 
torque the spacecraft about its axes for attitude control. The basic orbital 
maneuvers scheduled to be performed by the SPS for ATS-6 are initial station 
acquisition, east-west stationkeeping at 94 degrees west and 35 degrees east 
longitude, a transfer from 94 degrees west to 35 degrees east longitude (India) 
and the return to 105 degrees west longitude, and north-soutli stationkeeping 
as needed. 
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Figure 1-1. Alternate Orbital Control Jets (OCJ) 
Bar Thruster Configuration 
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Figure 1-2, Ion Engine Thmst Vector Geometry 
Relative to So^cecraft Axes 
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The continuous mode of thruster operations is used for all orbit control func- 
tions of the spacecraft. In this mode, tliruster ON times range from 1, 0 sec- 
ond to approximately 6 hours. The OFF time durations generally range from 
hours to days or months, so that each start is considered to occur at the am- 
bient hardware temperatur ^ . The other mode of operation is the pulse mode. 
This mode is used almost exclusively for attitude control functions on the 
spacecraft. It is generally characterized by very short ON times (.10 second 
(100 milliseconds) to 120 seconds), with relatively long OFF times (10 seconds 
to days). When orbital maneuvers are being performed, the Attitude Control 
System (ACS) stabilizes the spacecraft in the jet only mode. Attitude is main- 
tained by the attitude thrusters, and momentum wheels are not used for re- 
action torquing. 

As stated earlier, die orbital thrusting produces torques. Hence, it is not 
unlikely for an attitude thruster to fire in the pulse mode during an orbital 
maneuver. Thus, if thruster 7 is being used for a normal east-west station- 
keeping maneuver, thruster 3 or ll (the pitch thruster) may be fired in pulsed 
bursts to correct errors in pitch- The perturbations of the orbital elements 
arising from the attitude reaction thruster during orbital maneuvers have been 
estimated to be less than 1 percent of the desired maneuver. By deliberately 
misaligning the orbital thrusters, the thrust produced by the attitude thrusters 
will aid in accomplishing the desired maneuver. 

Figure 1-3 contains the nominal SPS performance data as a fiinction of tanl?: 
pressure. The data were reduced from tests conducted during the qualification 
and acceptance test phase. Appendix A contains a more detailed analysis of 
SPS performance characteristics. 
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1.2.2 Orbital Maneuvers With the Ion Engines 

Tile ATS-6 spiicecraft carries two identical Cesium bombardment ion engines 
for experimental use in maldng inclination and node corrections. Upon com- 
mand, the engines can deflect the thrust vector direction ± 7 degrees from the 
nominal direction. Thus, the ion engines can be positioned to thrust through 
the center of gravity. It is planned that the duty cycle during the ion engine 
experiment will be 7 hours of continuous operation for each 12-hour interval 
centered at the nodal crossings, covering a 1-month period. The thrust level 
of the ion engines is on the order of 0, 001 pound force, and the specific im- 
pulse is approximately 2500 seconds. 

1. 3 FUNCTIONAL SPECIFICATIONS FOR CNTRLF 

The original functional speciSeations for CNTRLF were written in June 1972. 
These specifications '‘';ere addressed In Reference 1, As a result of the de- 
cision to use SMSMAN as a base program instead of CONTRL, new functional 
specifications were written. In this section the new specifications are given 
as a reference to the reader. 

CNTRLF will compute maneuver control data for orbital and stationkeeping 
maneuvers to be performed by the ATS-6 spacecraft. It will be an offline 
program designed for use on the IBM 360 computers at the Goddard Space 
Flight Center (GSFC) and will be written in the FORTRAN IV language. 

Maneuver control data will be computed for the following SPS maneuvers: 
initial orbit correction, east-west stationkeeping (at 94 degrees west, 35 de- 
grees east longitude), the transfer to India and the return, and north-south 
stationkeeping. Maneuver control data will be computed also for the norih- 
south stationlceepirg maneuvers for the ion engine experiment. 

It will be possible to compute maneuvers in any one of four possible attitude 
confi^rations — reference, 3 -axis, spin and Sun, 
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1. Reference — This will be the normal operational attitude mode for 
the spacecraft when orbital maneuvers are being performed. Input 
attitude data will not be required. The program will compute the 
current attitude during the maneuver. 

2. 3-Axis — The three Euler angle rotations (yaw, roll, and pitch; will 
be entered and will remain constant throu^out the computations. 

3. Spin — This option will be used during initial acquisition (launch 
mode) when the spacecraft is spinning about its X-axis. The right 
ascension and declination of the spin axis will be the input data used 
to describe attitude. The attitude in inertial coordinates will be 
held fixed throughout the computations. Due to thruster orientation, 
only east-west maneuvers will be computed in this attitude mode. 

4. Sun — This option, like the spin option, will be used only during 
initial acquisition when the spacecraft' s X-axis is known to be 
aligned with and spinning about the Sun line. The input to CNTRLF 
will consist of the time of this alignment. The Sun declination and 
right ascension will be known a priori. The orientation of the 
ATS-6 spin axis wUi be set to the values at Sun sighting. 

The following models will be internal to CNTRLF; 

1. Hydrazine engine propulsion model 

2. Ion engine thrust perturbation model 

3. Finite burn integration model (fourth-order Adams-Moulton) 

The CNTRLF program will be able to operate using the following sources of 
orbit data: 

1. Tape or disk input (ORBl, EPHEM, GTDS sequential orbit file) 

2. Disk input (GTDS ORBIT file) 

3. Internal orbit generator (BPACKT/BPACKC) 
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Card input to the program will contain the following information; 


1. Initial spacecraft attitude 

2. Initial spacecraft orbit source and data, if necessary 

3. Desired maneuver 

4. Desired time to perform maneuver 

5. Engine to be used 

6. Engine and fuel tank initial states 
Y. Current engine calibration factors 

Program output (line printer) will contain the following information: 

1. Input card images 

2. Input summary 

3 . Propulsion system table 

4. Orbit and drift table 

5. Maneuver summary 

6. Control sheet 

7. Appropriate error messs^es 

The program will malffi the following nonrestrictive assumptions; 

1. Operating configuration parameters will not be specified (e. g. , 
heater turn-on times, latching value settings, coding of commands). 

2. The fli^t system data will be averaged into a homogeneous thruster 
model. This single model will contain thrust and specific impulse 
(ISP) data as a function of pressure, time, and temperature. A 
calibration factor will be used to adjust the homogeneous model 

for a specific thruster. 

3 . The internal program models will neglect the effect of attitude 
thrusting during orbital maneuvers. The effect will be included 
by proper adjustment of the thruster calibration factors. 
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The reference time will be Greenwich Mean Time (GMT) . 

Sngine cool down will not be modeled; thus, every bum will 
proceed as if the engine started at the ambient temperature. 



SECTION 2 - 3VIATHEMATICAL MODELS 


The mathematical models included in the CNTRLF program are described in 
this section according to the following categories: 

1. Orbital element systems 

2. Coordinate systems 

3. Maneuver options 

4. Hydrazine engine propulsion model 

5. Finite burn integration model (fourth-order Adams -Moulton) 

6. Ion englie thrust perturbation model 

7. Orbit sources before the burn 

No attempt is made to show how the models are integrated to obtain a working 
program. Reference 3 contains the necessary information on program struc- 
ture and interfaces. 

2.1 ORBITAL ELEMENTS 

The spacecraft orbit is specified by one of the following sets of generalized 
coordinates: 

1. Inertial position and velocity vectors 

2. Osculating Keplerian elements 

3. Injection or flight parameters 

4. Mean osculating elements 

2.1.1 Inertial Position and Velocity Vectors 

The inertial position and velocity vectors are defined as follows: 

R - (X , Y , Z) = the position or radius vector of the spacecraft from 
the center of the Earth expressed in the inertial coordinate system 

(X , Y , Z) = the velocity vector of the spacecraft expressed in the 
inertial coordinate system 
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The usual units of R and V components are Earth radial units (ERU) and 
ERU per second, respectively. One ERU equals the average radius of the 
Earth, or 6378,140 kilometers. 

2.1.2 Osculating Keplerian Elements 

The osculating Keplerian elements are defined as follows: 

a = semimajor axis 
e == eccentricity 
i = inclination 
m - mean anomaly 

argument of perigee 

S2 - right ascension of the ascending node 

The usual set of units for osculating Keplerian elements are either a in Idlo 
meters, e dimensionless, and i , m , tc , in degrees; or a in ERU, e 
dimensionless, and i , m , W , in radians. Figure 2-1 illustrates four 
of these elements, 

2.1.3 Injection or Flight Parameters 

The injection or flight parameters are defined as follows: 

0 = geodetic latitude of spacecraft subsateUite point in degrees 
{measured positive for north latitude) 

X = longitude of spacecraft subsatellite point in degrees (measured 
positive for longitude west of Greenwich) 

r = magnitude of the spacecraft radius vector in nautical miles 

V = magnihide of the spacecraft velocity vector in feet per second 

y - flight path elevation angle; the elevation of the inertial velocity 
vector in degrees (measured positive outward from the plane per- 
pendicular to the radius vector) 
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Figure 2-1 


OsculatiBg Keplerian Elements 







GL = flight path azimuth angle; tha azimuth of the inertial velocity vector 
projection onto the plane perpendicular to the radius vector in degrees 
(measured p:)sitive clockwise from east) 

Figure 2-2 illustrates these parameters. 

2.1.4 Mean Osculating Elements 

The mean osculating elements are used in Brouwer theory as an intermediate 
step in propagating the orbit. For details concerning mean elements, see 
Reference 5, 

2.1.5 Orbital Element Conversion 

Standard conversion formulas exist to convert elements of one system to those 
of another (Reference 6). These formulas have been programmed into a stand- 
ard set of software (Reference 5) which is used in many of the operational pro- 
grams at GSFC. Both CNTRLF and SMSMAN have combined these routines 
into a single package. This reduces the conversion task to a single call to a 
main driver. Subroutine EBMNT is used as a main driver of the standard 
conversion subroutines. The existing and the desired element types are each 
specified through ELEMNTs argument list by one of the following alphameric 
names; INERTIAL for inertial position and velocity vectors; OSCULATE for 
osculating Keplerian elements; INJECTION for injection parameters; and 
MEAN for Brouwer mean osculating elements. Thus, if the input parameter 
is INERTIAL and the output name is OSCULATE, ELEMNT calls those stand- 
ard routines which will convert inertial position and velocity vectors to oscu- 
lating Keplerian elements. 
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2*2 COORDINATE SYSTEMS 


The Earth-centered inertial (ECI) reference frame is defined by the following 
set of unit vectors ; 


= direction of mean equinox of date 

i xt 

1X1 

k. = direction of north celestial pole 
1 

When the SPS or ion engines are to be used for orbital maneuvers, the direc- 
tion of the AV vector or tlie thrust vector must be computed in tlie ECI frame. 
The unit vector obtained will be changing slowly as the spacecraft maintains 
a 3-axis stabilized condition, pointing at a location on the Earth. 


The calculation of the thrust vector in die ECI frame is done in a two-step 
process. First, the direction of thrust is transfr . from the spacecraft 
body axis frame to the local vertical (LV) frame (see Section 2.2. 1). The 
second transformation is from the LV frame to the ECI frame (see Sec- 
tion 2.2.2). 


The spacecraft body axis frame is defined by the following unit vectors (see 
Figure 2-3): 


-CV 

= parallel to mounting axis of solar array panel and positively directed 
away from optical axis of Polaris sensor 




= through center of parabolic reflector and center of Earth- viewing 
module; assumed positive in that direction 


The LV reference frame is defined by the following unit vectors (see Fig- 
ure 2-3): 


i^^ = east, parallel with the ecpiatorial plane 

^LV ~^V ^^LV * south 

= down, throu^ the Earth's center of mass 
JbV 
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Figure 2-3. Euler Eotation From Body Frame to LV Frame 




Table 2-1 gives the body coordinates for the SPS thrusters and ion engines 
aboard ATS-6. Note tiiat the east-west thrusters are canted, thus introducing 
a component of thrust away from the Earth for both eastward and westward 
thrusting. The cant angle has a magnitude of fi.575 degrees. The ion engines 
are nr oiuited on a pivot and are movable by command from the ground. The 
angle between the spacecraft’s Z-axis and the direction of thrust of the ion 
engine (C^) has been included to model this condition. 

The transformation from the spacecraft body axis frame to the ECI frame de- 
pends on the spacecraft attitude, the spacecraft orbital position, and the di- 
rection of the thrust relative to the spacecraft embedded axis. The ATS-6 
spacecraft attitude is entered into CNTRLF in one of the following modes 
(Section 2.2.3): 

1. Reference 

2, 3-Axis 

3. Spin 

4, Sun 

During the lifetime of ATS-6, it is expected that the above four methods of 
specifying attitude will cover all possible configurations in which an orbital 
maneuver could be performed. The user-specified orbit source (see Sec- 
tion 2,7) defines the initial orbit position. Table 2-1 defines the orientation 
of the thrust vector in the spacecraft body axis. Sections 2.2,1 and 2, 2. 2 de- 
scribe the transformation from the spacecraft body axis frame to the ECI frame 
assuming that the spacecraft attitude, orbital position, and embedded axes have 
been defined previously. Section 2,2.3 describes the four possible methods 
of specifying attitude Information, 
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Table 2-1. Body Coordinates for the Orbital Thrusters 
Onboard AT6-6 


HYDRAZINE ENGINE 

AV AND THRUST 
DIRECTION 

THRUSTER 

NUMBER 

THRUST VECTOR 
COMPONENTS IN 
BODY AXIS FRAME 


NOMINALLY WEST 

7 OR a 

II 

NOMINALLY EAST 

160R 16 

II 

NOMINALLY SOUTH 

BAND 13 

■Tb = 

NOMINALLY NORTH 

6 AND 14 


ION ENGINE 

SOUTH OR NORTH* 

NOT APPLICABLE 

<1 


•DEPENDS ON SIGN OF C 

a 

IS THE EAST-WEST THRUSTERS CANT ANGLE (A CONSTANT VALUE OF 5.575 DEGREESI 
C„ ISTHE ANGLE THE ION ENGINE MAKES WITH THE SPACECRAFTS Z-AXIS (NOMINALLY ±A5 DEGREES) 

ot — 
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2.2.1 Transformation From Spacecraft Body Axis Frame to LV Frame 

The spacecraft attitude is defined by tlie following set of three Euler rotations 
from th«f spacecraft body system to the LV frame; 

1. Yaw angle, 0 

2. Roll angle, 0 

3. Pitch angle, 6 

Figure 2-3 defines each of these individual rotations. 

Given the attitude (0 s 0 » 0 ) , the Euler transformation matrix used in the 
transformation from spacecraft body axis to the LV frame is defined as fol- 
lows: 


M- M, 

0 0 0 


where M 


0 


is the yaw angle transformation matrix 



cos 0 
sin 0 
0 


-sin 0 
cos 0 
0 


0 

0 

1 



is the roll angle transformation matrix 



1 

0 

0 


0 0 

cos 0 -sin 0 

sin 0 COB 0 


(2-1) 


(2-2) 
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and M. 
u 


is the pitch angle transformation matrix 



cos 0 

0 

sin 

0 

Mg = 

0 

1 

0 


-sin 0 

0 

cos 

0 


(2^3) 


Solving for M in the pitch, roE, yaw sequence, 


M = 


cos it cos 0 - sin (i sin if> sin 0 
sin li cos 0 + cos ili sin 0 sin 0 
-cos 0 Bin 0 


-sin ^ cos 0 cos ^ sin 6 + sin i/i sin 0 cos 0^ 

COB ij) COB 0 sin 0 sin 6 - cos 0 sin 0 cos 0 

sin 0 COB 0 cos 0 


(2-4) 


Defining T to be the unit vector in the direction of thrust in the LV frame 
^ & 

and T, the unit vector in the directiGn of thrust in the body axis frame, it 
b 

follows that 


T =M.T, (2-5) 

i b ' 

2.2.2 ’ Transformation From LV Frame to ECI Frame 

The azimuth and the elevation of the spacecraft position in the ECI frame are 
T). and ?7 , respectively (Figure 2-4). Defining T. to be the unit vector in 
the direction of thrust in tlie LV frame and T. to be the unit vector in the di- 
rection of thrust in the ECI frame, it follows that 

T =M T„ (2-6) 

i rj ^ * 


2“U 




Figure 2-4. Position of Satellite in Inertial 
Coordinate System 
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where 


M ^ 

n 


-sin 17 cos 7} sin 17 
X X A 


•cos 77 ^ cos 77^ 


cos ? 7 - sin ?7 sin 17 „ -sin rj cos 

i ±6 X a 


“COS ?7_ 


-Sin 77, 


(2-7) 


2.2.3 Modes of Attitude Operation 

In the normal mode of orbital thruster operation, the spacecraft attitude is 
held in the reference mode (Section 2.2, 3. 1). To model this mode, the pro- 
gram must calculate the three Euler angles (yaw, roll, and pitch) expected 
during the thrusting period. VS^en the satellite is pointing at a selected sta- 
tion on the Earth, the spacecraft will not be maneuvered Into the reference 
mode. In this mode of operation the three Euler an^es will be known a priori, 
and the angles will be entered as input, using the 3-axis program option (Sec- 
tion 2. 2, 3. 2). During initial acquisition when the spacecraft is spinning about 
its X-axis and is not 3-axis stabilized, the right ascension and declination of 
the spin axis will be loioivn directly or could be computed a priori from the 
time of aligmn^t with the Sun line. If orbital thrusting is requested during 
this period, CNTRLF will be nm in the spin or Sun modes of attitude opera- 
tions (Sections 2. 2. 3. 3 and 2. 3. 3.4, respectively). 

2 . 2 . 3 , 1 Reference Mode of Operation 

This is the normal mode in which attitude is maintained during orbital maneu- 
vers. The Euler angles of rv ’I and pitch are set to zero while the Eiiler angle 
of yaw is held so that the thrust vector will be, as closely as possible, either 
perpendicular or parallel to tlie orbit velocity vector. If external torques due 
to the orbital thrusting cause tlie spacecraft to deviate from this attiinide, the 
appropriate SPS thrusters are fired in the pulsed mode to reduce the error. 
The Euler angle of yaw is calculated so as to align the spacecraft X-axis as 



closely as possible witb the orbit velocity vector. Figure 2-5 defines the 
calculated yaw magnitude. The mathematical models used by the CNTKLF 
program to compute yaw are: 


E 


-cos B 
sin B 
0 


(2-B) 


where B = arctpn [R /R ] 

X y 

B "X component of spacecraft position vector In ECI frame 
R = y component of spacecraft position vector in ECI frame 


N = 


RX V 
|rx V 


/X 

T 


N X R 


ij) = arccos (T • E) 


(2-9) 


where V = satellite velocity vector in the ECI frame 
N - vector normal to orbit plane in ECI frame 
^ - vector tangential to ox-bit plane in ECI frame 

The three Euler an^es are Imown (yaw = calculated value, roll == 0, pitch - 0) 
as is the direction of the thrust in the body axis (see Table 2-1). Equation (2-5) 
is used to calculate the direction of thiust in the LV frame. I^nowing the 
spacecraft position in inertial space. Equation (2-6) is then used to calculate 
the direction of thrust in the ECI frame, 

2 . 2. 3 . 2 3 -Axis Mode of Operation 

At times the spacecraft will be tracking an Earth station or low-altitude space- 
craft during a time period in which an orbital thrusting mane r is required. 
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Figure 2-5* Definition of Parameters Used 
in the Heference Mode 


If the three Euler rotations at the desired time of the maneuver are Icnown, 
they can be entered directly into CNTBLP, The three Euler angles are known 
(yaw = input, roll = input, pitch - input) as is the direction of the thrust in the 
body axis (see Table 2-1). Equation (2-5) is used to calculate the direction 
of thrust in the LV frame. Knowing die spacecraft position in inertial space. 
Equation (2-6) is then used to calculate the direction of thrust in the ECI 
frame. 

2 . 2 . 3 . 3 Spin Stabilization Mode of Operation 

hi this mode, the spacecraft Is spinning slowly about the axis in which the east- 
west thrusters (7, 8, 15, and 16) lie (x-axis) anu the right ascension and 
declination of this axis in inertial space is Imowu. Therefore, the direction 
of thrust in the ECI frame can be solved for directly: 




T 

1 


cos Dec cos RA. 
cos Dec sin RA 
1 sin Dec J 


(2-10) 


/N 

vhere T. = the unit vector of thrust for thrusters 15 and 16 
^ (note that T. = -T. for thrusters 7 and 8) 

Dec = the user-specified spin axis declination 

RA = the user-specified spin axis right ascension 

2 . 2 . 3 . 4 Sun Stabilization Mode of Operation 

In this mode, the spacecraft is Imown to be aligned with and spinning about the 
Sun line. The time at which the spacecraft was aligned with the Sun line is 
entered into the program. A standard ATS routine SUNEPH (see Reference 5) 
is used to obtain the right ascension and declination of the Sun in inertial space 
at that time. By using Equation (2-10), the thrust direction in inertial space 
is computed. 
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2.3 MANEUVER OPTIONS 

IQ CNTRLF, thrusting requirements are computed for two general types of 
maneuvers (north-soth and east-west maneuvers) . A fiirther brealtdown can 
be made into classes which allow the user to model all orbital maneuvers ex- 
pected in the lifetime of ATS- 6. These classes are; 

• North- soutli maneuvers using the hydrazine thrusters (in pairs) and 
the ion engine, i. e, , inclination change maneuvers, node change 
maneuvers, and combined inclinaidon and node change maneuvers. 

• East- west maneuvers using only the hydrazine thrusters, i.e. , 
drift rate change maneuvers (drift rate goal specified by user) and 
stationkeeping maneuvers (drift rate goal computed by CNTRLF). 

CNTRLF models both nordi-south and east-west maneuvers in the following 
manner; 

1. Initial estimates of burn time and ignition time are developed 

2. The maneuver is modeled step-by-step until the estimated burn 
time is achieved. 

3. It is determined if the desired end conditions have been reached. 

If they have not, new estimates of burn time and ignition time are 
developed based on the error. Then, step 2 is repeated. 

This section is concerned with the determination of the initial estimates of bum 
time and ignition time. Those models used in the step-by-step iteration proc- 
ess are discussed in Sections 2.4 through 2. 6. 

2,3.1 Derivation of AV Magnitude for Drift Rate Change Maneuvers 

The drift rate goal and die ignition time have been either specified as input or 
derived as in Section 2. 3. 3. In either case, the drift rate to be achieved and 
die current drift rate are loiown (see Section 2. 3. 1. 1). A change in the drift 
rate is accomplished by adjusting the semimajor axis of the orbit. A positive 
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drift rate change (by definition, westward) is achieved by increasing the semi 
major axis; a negative drift change (by definition, eastward) is achieved by 
decreasing the semimajor axis. 

The change in the semimajor axis is accomplished by firing one of the east« 
west thrusters (7, 8, 15, or 16) in the orbit plane. The thrust is divided into 
a tangential component and a radial component (IJ^gure 2-6). hi general, the 
thrust vector causes a change in both the spacecraft's tangential and radial 
velocity vectors. Thus, the total change in the velocity vector (AV) due to 
thrusting in the orbit plane is defined as follows; 


AV = 



(2-Xl) 


where AV_ = AV component in orbital radial direction 
ix 

AV^ = AV component in orbital tangential direction 
AV = magnitude of AV vector 


AV, 

O^V-^ctan^ 


R 


(2-12) 


where = angle between AV_ and the direction of AV 
AV R 

The change in the semimajor axis due to thrusting can be represented by tlie 
following Lagrange Planetary equation (Reference 1); 




. 2e sin f . „ . 2a 

Aa AV_ H 

Z 2 R Tjr 

7?Vl “ e ' 


AV 


T 


(2-13) 
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CENTER OF EARTH 
ALONG RADIUS VECTOR) 


R = RADIAL COMPONENT 
T >= TANGENTIAL COMPONENT 


Figure 2-6, Thrust Vector Geometry in East-West 
Orbital Maneuvers 
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1 

where Aa = change in seminiajor axis (EHU) 

f = true anomaly (radians) 

a = semimajor axis (ERU) 

e = eccentriciiy (dimensionless) 

7} = mean motion = /a^ (radians per second) 

e 

AV_ = thruster velocity component (AV) in orbital radial direction 
K 

AV^ == thruster velocity component (AV) in orbital tangential direction 
r = position magnitude (ERU) 

3 2 

- Earth’s gravitational constant (62630.01853 ERU /sec ) 

The definition for r in elliptic orbits (see Reference 6) is given as 


a(l e^) 

1 + e cos f 


(2-14) 


The change in the semimajor axis Aa is also obtained by subtracting the cur- 
rent value of the semimajor axis from the desired value. Both values are 
calculated from the current and desired drift rates 


[86164. 1064 y]r1 


2/3 


a = 


2ir 1 


360/ J 


(2-15) 


where a - semimajor axis (ERU) 

(1 = Earth's gravitational constant (62630. 01853 ERU^/sec^) 
0 

X - drift rate — either current or desired (degree per day) 
where all of the parameters are defined in Equation (2-13). 


^One ERU = 6378.140 kilometers 
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Equation (2-13) can be simplified as follows: 


Aa “ 



2 

e 


f(e sin f) (1 + e cos f) AV^ 


(2-16) 


If we define 


C 


1 



2 

e 


e sin f 


(2-17) 


and 

C = ^ - (1 + e cos f) (2-18) 


Equation (2-16) can be rewritten as 


Aa-C^ AVj^ + Cg AV^ {2-19} 


Then, by definition, (see Figure 2-6) 


AV^ = AV cos 
AV^ = AV sin 


( 2 - 20 ) 


Substitutihg Equation (2-20) into Equation (2-19) and solving for AV produces 


AV = 


Aa 

Cl cos 0^^ 4- Cg sin 0^^ 


( 2 - 21 ) 
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2. 3. 1 . 1 Drift Rate Algorithm 

The accuracy to which the first estimate of AV can be calculated depends on 
how accurat^y the initial drift rate is Imown, In CNTRLF, the process for 
finding the initial drift is the same as that used in the acquisition table genera 
ation program for ATS-6, FACTBL (see Reference 7). In FACTBL, 49 sub- 
satellite points, equally spaced in time, are averaged to calculate the average 
subsatellite longitude over a sidereal day. The drift rate is then calculated by 
subtracting the current average longitude from the average longitude one side- 
real day later. The drift rate and the average longitude are computed in 
CNTRLF as follows: 

X. = Xav - Xav (2-23) 

1 1 


where X. = initial drift rate in degrees per sidereal day (westward 

^ drift is positive) 

X av^ , X = average subsatellite longitudes over a sidereal day at 
current time and one sidereal day later, respectively 
(westward longitudes are positive) 


Xav. 


48 

2 \+(X795.0855n) 


49 


(2-23) 


where X, is the initial subsatellite longitude and X. rmrr- v is the sub- 

1 ® i+(1795, 0855n) 

satellite longitude at the initial time plus n times 1795. Q855 seconds. Both X. 


and X. 


are calculated by converting the elements from the user- 


i+(1795.0855n) 

specified orbit source to Injection elements (degrees) (see Section 2. 1. 3). 


Once thrusting has started, the position and velocity vectors on the user- 
specified orbit source are no longer valid. Therefore, the method just de- 
scribed of calculating the drift rate is also no longer valid. The sole orbit 
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information available after tlmisting is that which comes from the finite bum 
integra.t’on model (see Section 2.5)* If this derived information is used to 
calculate iu Equation (2-23) and the orbit is propagated using the internal 
generator BPACKT (see Reference 11), an expression for drift rate throughout 
the thrusting period is formed as follows: 

X= (X. -X..) (2-24) 

Tin 

where X == drift rate after thrusting has started (degrees per sidereal day) 

■ 

X^ = drift rate obtained by propagating the Keplerian osculating ele- 
ments from tlie finite bum integrator using the rntemal generator 
BPACKT (degrees per sidereal day) 

X, = base drift rate as calculated in Equation (2-22) (degrees per 
^ sidereal day) 

# 

X.. = drift rate obtained by propagating the Keplerian osculating ele- 
ments from the user-specified orbit source at ignition time using 
the internal generator BPACKT (degrees per sidereal day) 

• ■ 

The term (X - X, . ) combines into one term the mismodeling of the internal 
1 11 

generator (BPACKT). Included in this term are all orbit perturbations not 

modeled in BPACKT (e. g. , the effects of the Sun and the Moon on the drift 

■ 

rate). Another way of interpreting Equation (2-24) is to assume that X. fixes 
the drift rate at ignition time and that any change in the semimajor axis will 
cause a corresponding change in the drift rate. X^ - X,. is the expression for 

the change in the drift rate resulting from a change in the semimajor axis. 

• * 

For small bums (^ 500 seconds), the mismodeling error (X^ - X^) will remain 
relatively constant. For burns of the hydrazine tlirusters aboard ATS-G longer 
than 500 seconds, the mismodeling error becomes relatively insignificant com- 
pared to the desired change in drift rate. If the effects of the Sun and the Moon 
were modeled in tlie internal generator BPACKT, the mismodeling error would 
approach zero and the two methods of computing the drift rate (before and 
during thrusting) would be equivalent. 
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2. 3. 1,2 Iteration Scheme for Drift Hate Maneuvers 


As stated earlier, CNTRLF is designed on the principle of obtaining an estimate 
of the burn time and then modeling the burn to verify that this estimate will pro- 
duce the desired maneuver goals. If the estimate does not produce the specified 
goals, an iteration scheme is applied to refine the estimate. For drift rate 
changes while thrusting in the reference attitude orientation, the maneuver 
can start at any time and the maneuver goal can always be achieved. In other 
words, although the maneuver may not be executed at the most fuel efficioat 
time (i. e. , at a apsis), the maneuver goal will be met. When tlumsting in the 
other possible attitude modes, however, there will be positions in the orbit 
where the maneuver goal cannot be met. No attempt is made in CNTHLF to 
refine the user-specified ignition time. Rather, the ignition time is held at 
its user-specified value, while the bum time is adjusted to meet the specified 
maneuver objectives. If the input attitude configuration is such that the thruster 
is not aligned with the velocity vector, a fuel optimum solution will not be 
found. In this case a warning mesisage is printed to alert the user, but the 
computation proceeds. The burn time is updated at each iteration step as 
follows: 

1. The drift rate algorithm (Section 2, S. 1. 1) is used to obtain a drift 
rate at the end of a burn. 

2. The semimajor axis corresponding to the drift rate after tlie bum 
is calculated using Equation (2-14). 

3. This semimajor axis is subtracted from the desired semimajor 
axis — also calculated by Equation (2-14) — to obtain an error term. 
The error term represents the magnitude of the error in the esti- 
mated burn Hme. 

4. Using Equation (2-21), this error in the semimajor axis is converted 
into an error in AV . 


2-24 



5 . 


The error in AV is then converted to an error in burn time by 
using the rocket equation; 


(AAV/(G ISP)) - 

AT = Wj — (2-25) 


where AT = the error time increment used to update the previous 
estimate of bum time (seconds) 

= the expected fuel weight at the end of the maneuver 
pounds 

AAV = the error in AV between the maneuver goal achieved 
and the desired maneuver goal (calculated in previous 
step) (feet per second) 

G = gravitational acceleration constant at Earth’s surface 
(32, 174 feet per second^) 

ISP = specific impulse of engine to be used (the es^ected 
ISP at the end of the burn) (seconds) 

W = rate of fuel usage for the engine to be used (pounds 
per second) i*e. , W = T/ISP where T equals the 
thrust of the engine to be used (pounds) (the expected 
thrust at the end of the burn) 

2.3.2 Derivation of AV Magnitude For Inclination/Node Change Maneuvers 

Inclination or node maneuvers (plane change maneuvers) take place at the in- 
tersection of the plane of the initial orbit and the plane of the final orbit (line 
of relative nodes). The initial orbit plane expressed in osculating Keplerian 
elements is obtained from the orbit source selected by the user. In general, 
the orbital vector read from the source will be at a time specified by the user. 
This time will not he the ignition time, as in the drift rate scheme (Sec- 
tion 2.3. 1.2), but will be an arbitrary time normally chosen very close to the 
desired maneuver time. The CNTRLF program will calculate an ignition 
time which will be as close to this arbitral^ time as possible (* one orbit 
period) and permit the north-south maneuver g lals to be met. The final orbit 
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plane is also eispressed in osculating Keplerian elements defined by the maneu- 
ver goals specified by the user. A change in the orbit plane due to thrusting 
can be represented by the following Lagrange equations (see Reference 1); 


dU r sin u 

dt " 2 r . 

na / 1 - e sm i 


(2-26) 


and 


di r cos u 

dt“T77=l^ 

na y 1 - e 


(2-27) 


where 


u = argument of latitude (argument of perigee plus true anom- 
aly) of maneuver point (radians) 

r - magnitude of orbit radius vector (ERU) 

a = semimajor a3ds (ERU) 

n - m 3an motion (radians per second) 

e = eccentricity (dimensionless) 

2 

W - orbit normal acceleration (ERU per second ) 
di, dJ2, dt = differentials of inclination, node, and time, respectively 


The Lagrange Equations (2-26) and (2-27) can be r nodified by replacing Wdt 
by dV (velocity differential) and by setting na' ^yT77r' equal to 
(orbital velocity). The modifications result in the following equations: 


dQ 


sin u 

: — rdv 

Vq sin 1 


(2-28) 


di = 


cos u 



dV 


(2-29) 
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Taking the partial derivatives of the changes in the node and the inclination 
with respect to the argument of latitude results in 


B(d^) _ cosu 
Vq sin i 


(2-30) 


au 


{2-31) 


Taking the partial derivatives of the changes in the node and the inclination 
with respect to the dV results in 


5(d^) _ sin u 
3(dV) “ Vg Bin i 


(2-32) 


3(di) _ cos u 


(2-33) 


By expanding the above results into a series, 


d(dfl) 


cos u dV 
sin i 


du + 


sin u 


V sin i 
0 


d(dV) 


(2-34) 


d(di) 


-sin u dV 



du + 


cos u 



d(dV) 


(2-35) 
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Solving Equations (2-34) and (2-35) simultaneously, using Cramer's rule re- 
iilts in 


du = [d(d S2) sin i cos u - d(di) sin u] (2-36) 

d(dV) = Vq [d(dQ) sin i sin u + d(di) cos u] (2-37) 

Approximate solutions of Equations (2-36) and (2-37) can be obtained if the 
finite differences are assumed to be equivalent to the differentials of these 
equations. The finite differences of Equations (2-36) and (2-37) take the 
forms 


Au = ” [EaA fisini^^ cosu]- (AAi sinu)} (2-38) 

AAV = V [ [ AA Osin i sin u] + (AAi cos u)} (2-39) 

u av 


where 


An = the finite change in the argument of latitude (radians) 

AAV= the finite change in the change in velocity increment (feet per 
second) 

AA0= the finite change in the change in node (radians) 

AAi = the finite change in the change in inclination (radians) 

= the average inclination (A i/2) (radians) 

= orbital velocity 
u = argument of latitude 
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2, 3. 2, 1 Calculation of the Time of Relative Nodes 

Solving Equations (2-28) and (2-29) simultaneously results in the following ex- 
pression for argument of latitude: 


u = arctan 


[ sin i dn 
( di j 


(2-40) 


where u = argument of latitude (argument of perigee plus true anomaly) 
(radians) 

i = inclination (radians) 
di, dO = differentials of inclination and node 


Talking finite differences, 


u = arctan 


( sin i 

av 

Ai 



(2-41) 


From the argument of latitude, the eccentric anomaly can be found 

f = u - CO (2-42) 

where f = true anomaly (radians) 

u = argument of latitude (radians) 

CO = argument of perigee (radians) 

Then, from Reference 6 


ECC - 2 arctan 


/X - e 
1 + e 


tan 



where ECC = eccentric anomaly (radians) 
e = eccentricity (dimensionless) 


(2-43) 
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The mean anomaly corresponding to tlie time of relative nodes is found from 
the eccentric anomaly. 


m = ECC - e sin (ECC) (2-44) 

where m = mean anomaly (radians) 

By subtracting this mean anomaly from the current mean anomaly and dividing 
by the mean motion of the orbit, tlie time ft'om the current time to the relative 
node is found. If the computed time is before the epoch of the orbit source, it 
is moved forward as many orbital periods necessary to place the time of rela- 
tive nodes after epoch. 

2.3.2. 2 Calculation of First Estimate of Burn Time and Ignition Time 

Knowing the time of relative nodes and an initial estimate of AV , the first 
estimates of the burn time and the ignition time can be calculated. The first 
estimate of the bum time is calculated using Equation (2-25). The values of 

9 

ISP and the fuel expenditure rate W are obtained at the start of the bum 
rather than at the end. Also the initial estimate of AV is used in place 
of AAV . The estimated bum time Is centered about the line of relative 
nodes to obtain a first estimate of the ignition time. 

2. 3. 2. 3 Integration Scheme for Inclination/Node Maneuvers 

An integration scheme is used with the inclination/node maneuvers to converge 
on a burn time (as with the drift rate change maneuvers) and an ignition time 
wMch will meet the user-specified maneuver goals. At each stage of file itera- 
tion, the burn time and the ignition time are updated according to the following 
procedure: 

1. An error is developed between the calculated final orbital inclination/ 
node and the desired inclination/node. 
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The errors in inclination and node become the AAi and AAQ » 
respectively, in Equations (2-38) and (2-39). 

3. Equation (2-39) is solved for AAV . The AAV is converted to 
an error in burn time using the Equation (2-25). 

4. Equation (2-38) is solved for Au . The Au is converted to an 
error in ignition time by dividing the results by the mean motion 
(n) of ihe orbit. 

2.3.3 StaHonlceeping 

Section 2, 3. 1 described the derivation of AV for drift rate change maneuvers. 
It was stated that the goal drift rate and the ignition time would be either user- 
specified or, m the case of east-west stationkeeping, calculated by CNTRLE. 
This section describes the calculation of the drift rate goal and the ignition time 
for east-west stationkeeping maneuvers. Section 2.3. 3. 1 discusses the station- 
keeping strategy incorporated in CNTRLF and Section 2. 3. 3,2 describes the 
specific targeting algorithm used. 

2 . 3 . 3 . 1 Stationke eping Strate gy 

The perturbing forces on an artificial statellite in orbit about the Earth are: 

1. The Earth ^s gravitational field 

2. The gravitational attractions of the Sun, the Moon, and the planets 

3. The Earth's atmosphere 

4. The Eai’th's magnetic field 

5. Solar radiation 

For circular, equatorial, synchronous orbits (well outside of the Earth's at- 
mosphere), the Earth's gravitational field and the gravitational attraction of 
the Sun and the Moon are the major perturbing forces. If it is assumed that 
the Earth is symmetrical about its polar axis, the magnitude of the Earth's po- 
tential field should be solely a function of the distance from the center of mass. 
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In general, iiie Earth's potential field departs slightly from that of a body 
having axial symmetry. The potential field varies not only as a function of 
distance from the Earth's center, but also as a function of longitude. The ac- 
celerated drift introduced in 24-hour orbits, due to the Earth’s anomalous 
potential, is described mathematically in Reference 9 as: 


V C, [F^ sinmX] + S, [-F, cos mX] (2-45) 
j^m Jtm Xm Am j \ / 


j2-m 


even 


where X = rate of change of drift rate (acceleration) due to the Earth 

not having axial symmetry (degrees per sidereal day^) 

X = subsatellite longitude (degrees east) 

C , S , = harmonic coefficients of the Earth's potential field 

j0m Am 

F „ = term to correct for inclination; 

Am 


F = — 
22 2 
a 


1 .xl2 

— (1 cos i) f 

2 i 


F 31 If cos i) + I sin^ i (1 + 3 cos i) 

1 2 a 


F =/Ml 


— (1 + cos i) 
2 


F = 

42 \ 4 

\a j 


1,^ .,2 7.2. 

— (1 + cos 1) - — sin 1 cos 1 (1 + cos 1) 


F. . = 


(4201 


44 I 4 


— (1 + cos i) 


a = semimajor axis (ERU) 
i = inclination (radians) 
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Included are the even A-m terms up to the higher order = 4 , m = 4 term* 

The constants C , and S * are measures of the amplitudes of the various 
j&m j&m 

spherical harmonic coefficients of the Earth’s potential field. These terms 
are referred to as tesseral harmonic coefficients. Values for these terms 
have been derived by examining the orbits of high altitude satellites. The fol" 
lowing values used by the Goddard Trajectory Determination System (GTDS) 
are also used in CNTHLF (Reference 10); 


1 

21 

C 

S 

2 

2 

.1586511 xio"^ 

8969932 X lo“ 

3 

1 

.2161875 XIO"^ 

o 

. 2571596 X lO" 

3 

3 

.1025055 XX0“® 

.1949036 XlO”^ 

4 

2 

.7739965 X lo""^ 

c 

.1515418 x-lo” 

4 

4 

-8 

-.3608512 X 10 

.6386659 X lo"^ 


Due to die tesseral acceleration, a synchronous satellite placed in an equa- 
tcrial orbit over a specified subsatellite point will drift slowly from its assigned 
slr.(ion. If left to drift, the satellite would drift toward a station on the Earth’s 
minor axis {see Figure 3-7). Like a pendulum, it would tlien reasonate about 
this axis. For this reason, active stationkeeping strategy is needed to main- 
tain the satellite position within some specified constraints. 

A scheme which has been used successfully with other operational synchronous 
satellites is to choose two boundaries, one on either side of a nominal station 
longitude, and to keep the spacecraft within these by exercising 

active thruster control. Figure 2-8 illustrates the sequence of events in nor- 
mal east-west stationkeeping. The satellite drifts slowly until it reaches the 
boundary that is between tiie normal station longitude and the closest minor 
axis. At this point, a maneuver must be performed, or the satellite will drift 
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outside of the boundary. The calculated maneuver gives the satellite the de- 
sired drift rate in the opposite direction, i.e. , toward the second hoimdary. 

The required impulse is calculated so that the effect of the tessera! acceler- 
ation reduces the drift rate to zero when the satellite reaches the second 
boundary. The continued action of the tes serai forces then starts the satellite 
back toward tlie maneuver boundary. When the first boundary has been again 
reached, another stationkeeping maneuver is performed, and the cycle is re- 
peated. Typically, this drift cycle (time between statiohkeeping maneuvers) 
is from 30 days to 80 days. In the scheme described above, calculation of 
the reqLiired impulse in CNTRLF is based on Equation (2-4=5) (perturbations 
due to tesseral acceleration only). However, as stated earlier, the gravita- 
tional attraction of the Sun and the Moon are also major perturbing forces. 
Because the impulse calculated by CNTRLF is correct only in the absence of 
a solar and a lunar gravitational field, the maneuvers calculated will result 
in minor overshoots or undershoots of the goal boundary. 

2 . 3 . 3 . 2 Tar geting Algorithm 

hi CNTRLF, a user can specify an east- west stationl^eeping maneuver in two 
ways, hi the first method, two average subsatellite longitude boundaries as 
well as a search start time, stop time, and time increment are input to CNTRLF. 
The program searches the specified orbit source until either of the two longi- 
tudes is crossed by the satellite or until the search stop time is reached. If 
a boundary crossing occurs, a maneuver is then oalculat<fd. If no boundary is 
crossed, an error message is printed and CNTRLF eidts. In the second 
method, the input consists of the time that the maneuver is to be performed 
and one longitude bountkiry. The input longitude boundary and the average 
longitude at the iiiput time define the boundaries in which the satellite is to be 
constrained (see Figure 2-8). The average longitude calculated at the input 
time serves as the maneuver boundary, while the input longitude defines the 
goal boundary. With the time of the maneuver now known (either the result of a 
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search or entered directly, a goal drift rate is determined., Based on the con- 
cepts developed in Section 2. 3. 3. 1, the following algorithm is used to calculate 
a goal drift rate 


1. Determine within which quadrant the maneuver falls 


2 . 


Quadrant 


Boundaries 


11.5 degrees west - 75 degrees east 
75 degrees east - 161. 8 degrees east 
161, 8 degrees east - 105. 5 degrees west 

105. 5 degrees west - 11. 5 degrees west 


IQUADCl) 

IQUAD(2) 

IQUAD(3) 

IQUAD(4) 

If within quadrant 1 or 3, the calculation will result in either a 
zero or a positive drift rate 


3. If within quadrant 2 or 4» the calculation will result in either a 
zero or a negative drift rate 

4. If the maneuver boundary is closer to the minor axis than the goal 
boundary, the drift rate goal is calculated as follows : 


X= ±y2XAX 


(2-46) 


where X = calculated drift rate (degrees per day); sign deter- 
mined by quadrant in which satellite lies 

X = tesseral acceleration calculated using Equation (2-45); 
longitude used in Equation (2-45) is that of the desired 
desired station— center of ihe tolerance box (degrees 
per day^) 

AX = distance between the two longitudinal constraints 
(degrees) 

If the maneuver boimdaiy is farther from the minor axis than the 
goal boundary, the drift rate goal is zero. 
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2.4 HYDRAZINE ENGINE PROPULSION MODEL 


The hydrazine engine propulsion model includes the models for the following 
procedures: 

1, For a time segment of At seconds (continuous firing mode), com- 
puting the weight of fuel esqaended and update the tanlc pressure 
and spacecraft weight as a result of diis fuel expenditure. 

2, For a time segment of At seconds, computing the effective im- 
pulse and AV delivered by die control jet 

3, Interpolating the control jet calibration tables as necessary 
2.4,1 Fuel Weight Expenditure 

The hydrazine propulsion model used in CNTRLi? computes the weight AW of 
fuel expended during At seconds of continuous firing, 

AW + W.) A t (2-47) 

where AW = the weight of fuel expended during At seconds of continuous 
firing (pounds) 

W.^^ = rate of fuel expenditure at the end of At time segment; obtained 
by interpolating the W calibration tables (pounds per second) 

W. = rate of fuel e35>enditure at the start of the At time segment; 
obtained by interpolating the W calibration tables (pounds per 
second) 

The problem of computing AW reduces to interpolating the W calibration 
tables. These tables (see Section 2.4.3) express W as a function of engine 
type, pressure, and time. The engine type and time are known. The pressure 
selected for the interpolation will be the average pressure P between the 
initial pressure at the start of the time segment and che final pressure 
P^ at the end of the time segment. 
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(2-48) 


p = 
av 


P + Pn 
0 f 


wliere P 

av 



= average pressure in time segment (pounds per square inoh ab- 
solute (PSIA.)) 

= initial pressure at start of the time segment (PSIA) 

= final pressure at the end of the time segment (PSIA) 


The pressure at the beginning of the time segment is Icnown. Using Boyle’s 
law, the pressure at the end of tiie time segment, P^ , is expressed as follows t 


PCON 

^ Wq* AW 


(2-49) 


where P^ ~ final pressure at the end of the time segment (PSIA) 

PCON = pres sure -volume constant of the fuel pressurant given by 


PCON = Pj 


^WFULL - W 


fuel. 


P = tank pressure at the start of the maneuver — not of the time 
segment (PSIA); input to CNTRLF program 

WFULL = weight of fuel required to fill the tank completely (pounds); 

data in spacecraft calibration data base (see Section 2. 4. 3) 

Wfaelj “ wei^t of fuel at start of the maneuver (pounds); input to 
CNTRLF program 

Wq* = relative pressurant volume at the start of a time segment 
(pounds) given by 


W^* = WFULL 


W 


fuel 


0 


(volume is proportional to weight) 

^fuelo " start of a time segment (pounds); Internal 

variable maintained by CNTRLF 

AW = weight of foel expended (pounds) 
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Equation (2-49) reflects tlie decrease in tanl?: pressure due to ti\e expanding 
volume of pressuiant as fuel is depleted from the tanl^. 

To calculate the fuel expended in a time segment, Equation (2-47) indicates that 
die rate of fuel expenditure curves must be interpolated. To interpolate, the 
average pressure during the time segment must be ’ 'ti. is a 

function of the final pressure which is a function of the weigh expended. In 
CNTRLF, AW is determined iteratively using Equations (2-47), (2-48), and 
(2-40) as follows: 

1. An initial estimate of AW is made: AW = 0. 0; thus, P^ = P. 

and P = P. . 
av 1 

2. Using Equation (2-47), AW is solved using P - P. to interpo- 

fl.V 1 

late the W calibration tables. The interpolated W value is cor- 
rected for individual jet performance by multiplying tlie value by 
a user-specified calibration factor. 

3. A comparison is made between the AW estimated and the AW 
calculated. If the difference is less than 0. 00001 pound, the cor- 
rect AW and P have been found, 

av 

4. If AW estimated and AW calculated differ by more than 

0. 00001 pound, use as a new estimate of AW the previously cal- 
culated value. This value of AW is used in Equation (2-49) to 
determine a new value for the final pressure at the end of tlie time 
step P^ . All values in Equation (2-49) are known at the start of 
the time step except AW . Now, solve Equation (2-48) for a new 
estimate of P 

av 

5. Next, Equation (2-47) is solved for AW using the F value cal- 
culated using Equation (2-48) to interpolate the W calibration 
tables. The interpolated W value is corrected for individual jet 
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performance by multiplying the value by a user-specified calibra- 
tion factor. 

6, Step 4 in which, a comparison is made between the estimated AW 
and tlie calculated AW is repeated. 

In the above iteration, both the fuel expended in a At time segment and the 
average pressure P during that time segment are calculated. If after 
20 iterations the estimated AW and computed AW differ by greater than 
0. 00001 pound, an error message is printed and the final calculated value of 
AW is used. 

2.4.2 Effective Impulse and AV 

Once the fuel expended and the average pressure during a At time segment 
have been determined, the following calculations are made: 

N cos Qi 

EFFIMP = 2 (T, + T-) At (2-50) 

ii 12 

where EFFIMP = effective impulse (pound/seconds) 

N = number of thrusters used; 

= 1, if an east-west static nkeeping iuaiitju/er is performed 
= 2, if a north-south maneuver is performed 

d = cant angle which introduces a radial component of thrust for 
both eastward and westward thrusting; for north-south ma- 
neuvers the cant angle is set to zero degree 

thmst at the beginning of the AT time increment (pounds) 

thrust at end of AT time increment (pounds]. 

time segment in which calculation is to be made; generally 
lies between 1 second and 100 seconds 

T and T are found by interpolating the thrust calibration table (see Sec- 
1 ^ 

tion 2. 4. 3) at the beginning and end of the At time increment and with pressure 
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equal to the average pressure durmg the tiipe step . The interpolated 
thrust is corrected for individual jet performance by multiplying the value by 
a user-specified calibration factor. 

AVEISP = EFFIMP/AW (2-51) 


where AVESP = avex’age ISP over time interval At 

EFFIMP - effective impulse from Equation (2-50) 

AW = fuel suspended during At time segment; calculated iterati- 
vely (see Section 2.4. 1) 


DELVE == AVEISP G 



(2-52) 


where DELVL = AV expended due to thrusting duriug At time segment (feet 
per second) 


G = acceleration of gravity at Earth's surface (32.174 feet per 
second^) 


- spacecraft's weight before incremeptal burn; the weight at 
the start of tlie burn is put into CNTRLF; the spacecraft's 
weight is recalculated after each At time segment (pounds) 


= spacecraft's weight after incremental burn (pounds); calcu- 
lated by the expression - AW ; used as W^ value in 

the next time segment 


2.4.3 Thruster Model Data Base 

The jet calibration tables are stored internally in CNTELF. The data is in 
tabular foim and consists of the best available thruster data from ground tests 
(see Appendix A). The following list describes the data required by the engine 
models of CNTRLF. 
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Syrnbol 

Definition 

Units 

G 

Gravitational constant on Earth’s surface 

Feet per secor 
squared 

0! 

Hydrazine engine cant angle 

Radians 

THRION 

Thrust of ion engine 

Pounds 

XISP 

ISP of ion engine 

Seconds 

N 

Number of hydrazine engines modeled; a 
maximum of two are allowed. CNTRLF 
presentLy uses one model with the individ- 
ual thrusters varying from tl'+e model by a 
user-specified calibration factor 

Not applicable 

NOTMPS 

Number of tanlt temperatures at which hy- 
drazine thruster data is available; FDATA 
uses at least three temperature curves 

Not applicable 

TABTMP 

Inlet temperatures at which thruster data 
is available 

Degrees 

Fahrenheit 

TAB VOL 

Array of weights of fuel Ln tanks when com- 
pletely full (WFULL). There is one value 
for each inlet temperature 

Pounds 

NOPHES^ 

Number of pressure data points on indi- 
vidual temperature cuives 

Not applicable 

TABPRS 

Pressures corresponding to points on 
temperature curves 

Pounds per 
square inch 

NOTIME ^ 

Number of time data points on individual 
temperature curves 

Not applicable 

TABTIM 

Time corresponding to points on tempera- 
ture curve 

Seconds 

TABBP 

A four-dimensional array representing 
specific impulse at a specified time, pres- 
sure, and inlet temperature for a specific 
engine model 

Seconds 


^The number of time points and pressure points must be the same for each 
inlet temperature curve. The values of time and pressure at tliese points 
can be different. 
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Symbol 

Definition 

Unite 

TABF 

A four-dimensional array representing thrust 
at a specified time, pressure, and inlet tem- 
perature for a specific engine model 

Pounds 

TSEG 

Maximum burn time of a hydrazine thruster 

Minutes 

DELSEG 

Minimum allowable time between successive 
actuations of the hydrazine thrusters 

Minutes 


As new data becomes available, an internal subroutine (FDATA) will be mod- 
ified and recompiled. Currently, ATS-6 flight test data are incorporated in 
the program. The current data and format are shown in Fi^re 2-9. 

Based on an input tank temperature, the data in TABVOL are Inteipolated (see 
Equation (2-54)) to obtain a current value for the amount of fuel required to 
completely fill the tanks (WFULL). Based on an input tank temperature, the 
data in TABISP and TABF are interpolated (see Equation (2-54)) to obtain data 
describing ISP and thrust at that temperature. The thrust data is divided by 
the ISP data on a point-by-point basis to obtain data on the rate of fuel esqiend- 
iture. 


W = 


T 

ISP 


(2-53) 


where W = rate of fuel expenditure at a specified time and pressure (pounds 
per second) 

T = thrust at the same time and pressure (pounds) 

ISP = ISP at the same time and pressure (seconds) 

All data points are indexed by type of jet mc^eied, time, and pressure. The 

* 

fuel expenditure rate data W and the thrust data T are used in the calculation 
of fuel expended AW and average pressure P^^ during an incremental time 
step At (see Section 2, 4, 1). Specific values for W and T at a given time 
and pressure are found by inteipolation, A linear interpolation scheme is used 
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Figure 2-9. Current ATS-6 Thruster Data at 
Three Temperatures 
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(see Equation (2-54)) unless the current p^ressure falls within the three lowest 
pressure points; if this happens, a quadratic interpolation method is used. 

The data are interpolated first witli respect to pressure and tiien time using 
the following equation; 


F. 


X 




(^xu " ^xl) 


(2-54) 


e- 




(2-55) 


where X = independent variable 

Xt = lower independent variable data point 

X^ = upper independent variable data point 

F = lower dependent variable data point 
aL 

F = upper dependent variable data point 
AU 

F = inteipolated answer 
X 

In the above scheme, the independent variables are pressure and time, while 
tlie fuel esqjendxture rate and the thrust are dependent variables. 

Included in tlie list of data required by the CNTRLF engine models are the ion 
engine parameters. These parameters are constant and do not vary as a func- 
tion of time and pressure as do tiie hydrazine engine parameters. 

2.5 FIJ^TE BURN INTEGRATION MODEL 

Once the amount of fuel expended and AV magnitude produced during an in- 
cremental burn are calculated, then the changes in the orbital elements due to 
the thrusting of a hydrazine engine are calculated. To calculate these changes 
requires Imowledge only of the inertial position and velocity vectors at the 
start of the bum segment and of the type of engine being used. The Earth is 
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treated, a point source and a standard numerical technique is used to in- 
tegrate the R equation, which is 



+ AAV T 


(2-56) 


where R„ = final acceleration vector (ERU per second ) 

I 

jLt,= gravitat’ nnal attraction constant of the Earth (32. 174 feet per 
second*'^ 

^ = calculated position vector representing the spacecraft position 
at the end of a calculation time step (ERU); at ignition time, this 
value is obtained from the user-specified orbit source 

AAV = rate of change of AV during the incremental time step (feet 
per second) 

T - thrust direction unit vector in inertial coordinates; calculated 
each AT time step using the transformation described in 
in Section 2,2 

The specific numerical technique used in the fourth-order A dams -Moulton 
predictor-corrector integration scheme. To start the predictor-corrector, a 
modified second-order Runge-Kutta method is used. In this method, the fol- 
lowing equations, which are integrals of the R equation, are solved; 


u R 

^ ave 

y = V. -- — I^AT + AVT 
f X ave 

it 


R, = R. + V AT + DVR T AT 
f X ave ave 


(2-57) 

(2-58) 


where 


j 


V„ = velocity vector of the spacecraft at the end of the time step 
(ERU per second) 

R^ = position vector of the spacecraft at the end of tlie time step (ERU) 
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^ 1 

V = estimated velocity vector in the middle of the time step (ERU 
per second) 

"R = estimated position vector^ in the middle of the time step (ERU) 

$ = estimated thrust direction vector in the middle of the time step 

(unit vector) 

*V, = initial velocity vector at the heginning of the time step (ERU per 
second) 

R. = intial position vector at the beginning of the time step (ERU) 

At = incremental time step (seconds) 

AV change in orbital velocity due to tiirusting calculated using 
Equation (2-52) (ERU per second) 

DVR = Integral of AV defined by the following equation (with all vari- 
ables defined as in Section 2.4«2) 


DVR = AVEISP G 



(2-59) 


The starter method is used only for three At time intervals before the Adams- 
Moulton method of numeiical integration is initiated. If the At time interval 
changes (e. g. , at the end of the maneuver), the Adams-Moulton method is re- 
started with the modified second-order Runge-Kutta, 


1 

The estimated velocity, position, and thrust directions are calculated by 
advancing the mean anomsily of the os culating elements at the end of the pre- 
vious time step by half a time step, then converting the resulting osculating 
elements into inertial elements. 
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The fourth-order Adams -Moulton predictor-corrector is defined mathemati- 
cally as follows: 

Adams-Moulton predictor: 


PV - V. +~ (56A. - 59A. , + 37A. . - 9A, -) 

1 24 X 1^1 1-2 1-3^ 

?R =?. (55V. - Se'V. , + 37V! „ - 9V. _) 

1 24 1 1-1 x-2 1-3" 

Adams-Moulton corrector; 


—3^ AT 

CR = R. + ^ (9PV + 19V. - 5V. + V. „) 

1 24 1 x-1 1-2 


(2-60) 


(2-61) 


where = predicted acceleration calculated using the R Equa- 

tion (2-56) (ERU per second^) 

PV, PR = predicted velocity and position vectors^ of the spacecraft 
at the end of the incremental time step (ERU per second) 

CV, ^ = corrected velocity and position vectors^ at the end, of the 
incremental time step (ERU per second) 

V., R. = acceleration, velocity, and position vectors of the space- 
craft at the beginning of the time step 


A. 

i-x 


I 




acceleration and velocity vectors of the spacecraft at 
either one, two, or three time steps previous to the be- 
ginning of the current time step 


The predicted R and calculated R should not differ by more than a toler- 
ance of 1 X 10“® ERU . If they do an iteration scheme is entered until they 
converge within this tolerance. 
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CNTRLF maintains both an acceleration and velocity array where the following 
substitutions talie place at the end of every time step: 


A. „ 


i-3 

i-2 


„ 

V. « 


1-3 

i— 2 



A. „ 

= A 

i-2 

i-1 




= V 

*i-2 

1-1 



A. , 

= A. 

i-1 

I 


- 


= V. 

1-1 

1 


The final position and velocity vectors are transformed into Keplerian osculat- 
ing elements using the method discussed in Section 2. 1. 5„ 

2. 6 ION ENGINE THRUST PERTURBATION MODEL 

Sections 2.4 and 2.5 describe the effects of thrusting the hydrazine engines for 
an incremental time step. The following section describes the effects of thrust- 
ing the Cesium bombardment engine for an incremental time step. 

The thrust level of the ion engines onboard ATS-6 is on the order of 0. 001 pound 
force, and the specific impulse is about twenty-five hundred seconds. The 
thrust does not vary as a function of time or pressure and is only slightly 
greater (by a factor of two or three) than the external forces acting to perturb 
the orbit. Due to this low nonvarying thrust, the hydrazine engine propulsion 
models described in Sections 2.4 and 2.5 are not valid for the ion engines. 
Instead, a method which computes the change in orbital elements caused by 
continuous low-level thrusting over an ejctended period time is required. 

Similar to the method used for the hydrazine engines this computational method 


2-50 



subdivides the maneuvers into incremental time segments of At duration. 
Computed velociiy and position increments are vectorially added to the posi- 
tion and velocity vectors read from the initial orbit source at each computa- 
tional time step. As a result, the firing of ion engines is treated as another 
perturbation of the original orbit. To compute the velocity and position in- 
crements, the following recursive formulas are used: 


AV =AV +KF 

X ji X 


AV =AV +KF 

y y y 


(2-62) 


AV =^AV +ICF 
z z z 


AX==AX+ AV --^KF At 
X 2 X 


AY=AY+ AV “irKF At 
y 2 y/ 


(2-63) 


AZ = AZ + [AV --^KF 1 At 
z 2 21 


where = velocity increments (ERU per second) 

AX, AY, AZ = position increments (ERU) 

F^, F , F^ = thrust components obtained by multiplying the imit 
thrust vector (as computed in Section 2.2) by the 


X y z 


constant ion engine thrust level 
K = value of At/M (seconds per slugs) 
At = incremental time step (seconds) 

M - mass of the spacecraft (slugs) 


M = W,-^KAt 


I ISP. 


(2-64) 


ION 
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where 


M = mass of the spa> -^craft (slugs) 

Wj = weight of spacecraft at start of maneuver (pounds) 
'^lON ~ engine (0, 001 pound) 

ISP „ = specific impulse of the ion engine (2500 seconds) 

lvJ(N 

X =s number of time steps since ignition time 
At = incremental time step (seconds) 

The equations are used recursively as follows: 

At t = ignition time + At 


AV = 0 

1 M IVi 

AK = 0 -i- (aV --i ^ f) At "i — At^ 
2 M 2 hi 


Substituting the equation F = Ma , one obtains 


AV, = aAt 
1 

AX, =iaAt^ 

1 M 


At t = ignition time + 2 At 

AX^ = AX^.(AV^-|fp)Ata^(2At)2 
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Substitutiiig tli 0 equation F = Ma and setting 2 At = At_^_ * - , one obtains 

^OX AXj 


AVg a ^^total 


AXg 


Applying the same analysis to any other time (ignition time + x At ), one obtains 

AV^ = a At^Q^^j^ 

^^^TOTAL 


where AV = incremental AV to be added to velocity vector at time x 

AX - incremental change in position to be added to position 
vector at time x 


a - acceleration due to tiirusting 
At^OTAL ” incremental time (xAt) 


At each time step, the equations reduce to the standard set of velocity and 
position equations in a constant acceleration field. Also at each time step, 
the calculated position and velocity increments are added vectorially to die 
position and velocity vectors from die specified orbit source. The inertial 
orbital elements (Cartesian) formed are then converted to osculating Keplerian 
elements as described in Section 2. 1.5. The converted osculating Keplerian 
elements form the burnout elements from thrusting widi the ion engines for 
AT seconds. Note that the required thrust and ISP parameters are stored in- 
ternally in CNTRLF (see Section 2.4.3), 
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2. 7 ORBIT SOURCES BEFORE THE BURN 

CNTRLF obtains its initial inertial position and velocity vectors (converted to 
osculating Kepleriaii elements) before a bum by reading the standard ephemeris 
sources available at GSFC. If these are not available, CNTRLF can propa- 
gate an ox'bit using the internal analytic generator BPACKT. The values read 
or propagated are used to determine the initial orbital elements (both osculating 
Kepleiian and inertial position and velocity vector) at ignition time. These ele- 
ments are used by the finite burn integrator (see Section 2. 5) and also by the 
drift rate algorithm in calculating the initial drift rate (see Section 2.3. 1.1), 
The initial orbit source is also used for all ion engine calculations (see Sec- 
tion 2.6). When the orbit source is external to CNTRLF, it reside? either on 
a tape or on a disk and is in the form of an EPHEM file, ORBl file, or GTDS 
ORBIT file. 

2.7.1 Internal Orbit Source 

CNTRLF contains an analytic propagator based on Brouwer theory (Refer- 
ence 10). The internal generator is subroutine BPACKT. This subroutine 
serves as an executive routine to the standard ATS subroutine BPACKC 
(Reference 5), The propagator models the zonal and tesseral spherical 
harmonic terms which are used to describe variations in tlie Earth’s gravi- 
tational field. Not included are the important Sun and Moon gravitational 
perturbations. If none of the standard orbit sources at GSFC is available, 
the internal generator can serve as the initial orbit source. In this mode the 
user enters known epoch elements at a given epoch time. The elements are 
either Keplerian osculating elements, injection elements, or inertial position 
and velocity vectors (see Section 2.1), The input elements are converted to 
osculating Keplerian elements and then to Brouwer mean elements. The 
Brouwer mean elements are used to propagate to spacecraft orbit to some 
arbitrary time after or before epoch. The mean elements are then converted 
bad? into osculating Keplerian elements. The primary purpose for the internal 
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generator is to calculate tiie drift rate after thrusting has started (see Sec- 
tion 2. 3* 1). It is used in the drift rate calculation regardless of the initial 

orbit source. 

2.7,3 Standard Orbit Sources at G5FC 

At GSFC there are two large scale orbit determination systems: GTDS (Ref- 
erence 11) and the Definitive Orbit Determination System (DODS) (Refer- 
ence 12). GTDS produces as output EPHEM files, ORBl files, and GTDS 
ORBIT files. DODS has only the EPHEM and ORBl capacity. The EPHEM 
and ORBl files are sequential data sets containing epoch data and orbital po- 
sition and velocity elements in the form of vector. Note tliat EPHEM and 
ORBl files contain essentially the same information with different units. The 
GTDS ORBIT file is a direct-access data set which can contain the trajectories 
of multiple orbits. The user specifies which orbit is of interest through a 
level parameter. The data records in the file are read and interpolated using 
standard subroutines available for this purpose (Reference 13). In general, 
after tracking ATS-6 for 24 hours, the range and range rate data obtained are 
reduced Into an orbit by the orbit computation engineer using either GTDS or 
DODS. Once an orbit has been determined, it is propagated over a month's 
span and the results are output as one of the files described above. Using the 
standard subroutines, CNTRLE reads the header record which contains epoch 
information and the appropriate data record which contains the position and 
velocity vector at ignition time. If the vectors at the desired time are not 
contained witiiin the file or if there is an error in reading the file, the epoch 
information read is then used to initialme the internal orbit generator. In 
this case, an error message is printed and the internal orbit generator is 
automatically used as the initial orbit source. 
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APPENDIX A - ATS-6 HYDRAZINE PROPULSION SYSTEM DATA 


A nece$sary step in tlie development of the ATS-6 orbital maneuver control 
program (CNTRDF) is the acquisition and incorporation in the program of a 
data set which accurately defines the thrust and specific impulse of the hydra- 
zine tlirusters as a function of propellant pressure and temperature. On 
October 12, 1972, the author visited Fairchild Industries and met with 
the propulsion personnel* Ail data available at that time including the proto- 
lype system test data was provided. These data have been analyzed according 
to the following criteria: 

• Establishing the most appropriate format for utilization of the 
flight system data when these data become available 

• Establishing ground rules for math models of the hydrazine system 
in the CNTRLF program 

The analysis is summarized below together with the results, conclusions, and 
recommendations made on tlie basis of these findings. 

A. 1 PURPOSE OF ANALYSIS 

The purpose of the analysis was to determine the most appropriate format for 
using the test data. This will establish ground rules for modeling the hydrazine 
thruster performance. Consideration was given to the following questions: 

1. Which test data should be used? 

a. Individual thruster performance data 

b. System test data 

c. Both types of data 


fKECEDffiG page BLdEK NOT 



2 . 


How should the data be used? 


a. Average aH thrusters using a calibratiou factor to discrimi- 
nate between thrusters 

b. Model each thruster individually 

3. What thrust and ]BP models are dictated by the data? 

a. Time and pressure dependency 

b. Temperature dependency 

c. A combination of the above 

A. 2 DATA USED 

The type of tests and the data used in each follow. 

1. Performance Mapping Tests — Each thruster was tested at a matrix 
of selected inlet temperatures and tank pressures. The thruster 
burn time was either 200 or 300 seconds for steady state thrusting. 
Necessary data were taken to allow thrust computation, 

2. Individual Thinister Tests — Each thruster was tested at a constant 
inlet temperature of 70°F and varying tank pressures. The burn 
time of the thruster was 200 seconds. Necessary data were taken 
to allow both ISP and thrust computation. 

3. ATS-6 Sustained Firing Tests — One time engine tests were made 

of thrust and ISP at 1-, 6-, and IS-hour operating times. The inlet 
temperature and tanic pressure were held constant. 

All tests were conducted so &at the calculation error in iSP and thrust would 
be less than ±5 percent. 
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A. 3 TESTING CONDITIONS 


Pressure and temperature was read from a strip chart or, if they varied ra- 
pidly, an oscillograph was used. Calibration factors corrected equipment 
biases. In steady state flow the thrust was calculated from the test data by 


T == P 

c 



C 


f 


where T = the thrust in pounds 

P = the measured chamber pressure. It is converted from PSIG to 
^ PSIA. An assumption is tliat tank pressure is a monotonicly in- 
creasing function of chamber pressure 

2 

A^^ = measured thrust area corrected for thermal expansion (inches ) 

Cj = thrust constant to be determined by the project office. One 
way of obtaining C^ is as follows: 


Another expression for thrust, T , is 


T=6a C^ 
tc 


where 6 = propellant density 
C = exit velocity 


Thus 


c tc f 


6a^ C^ 
tc 


which reduces to 


P C 
c 


f 


6 
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Imowing Bernoulli's Equation 



it can be concluded C„ 2 . 

I 

In tests like the individual tliruster tesis, propellant flow W'as measured by 
sight glass allowing the following equations to be solved; 

W " W 


where ISP = the specific impulse in seconds 
T = the thrust in pounds 

W = the rate of weight flow in pounds per second 

w=-~ 

At 


where Ah = the sight glass deflection (centimeters) 

3 

K = the sight glass constant (in per centimeter) 
6 = the propellant density 
At = the time between readings of the sight glass 



A. 4 ANALYSIS PERFORMED 




To perforin the analysis, the data were graphed with tolerance lines and aver- 
age data lines drawn, when appropriate. The following illustrations were gen- 
erated from ihese graphs to compare thruster performance: 

Table A-1. Repeatability 

Table A-2. Individual Thruster Modeling 

Table A-3. Average Thruster Modeling 

Table A-4. Comparison of Individual Thruster Modeling to Average 
Thruster Modeling 

Table A-5, Error in Burn Time Due to Worst Cases 

Table A- 6. Final Calibration values 

Figure A-1, Final Thruster Model at 70°F 

Explanatory notes are presented with each of these illustrations. 

A. 5 CONCLUSIONS 

The following conclusions were drawn from the tables. 

1. Use ail data. Both tlie individual thruster test data and the per- 
formance mapping data should be used in the thrusting model due 
to the good repeatability of die data. 

2. Model a single thruster. By averaging all the data, a single 
thruster model can be developed. Tliis model is no worse thsui the 
worst case when an individual thruster is modeled. The model 
will give from 3 to 4 percent errors in burn time. 

3. Use calibration factors witii the model. The individual modeling 
shows on the average a 1 percent improvement over the average 
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modeling. Calibration factors will result in a similar improve- 
ment. The calibration factor will take the following form: 




T 

ave 


where T = the thrust for thruster X 

-?C 

= the calibration factor /1 00 

"^ave ~ average value of thrust for a nomogeneous model 

4. Model thrust and ISP as a function of times pressure and inlet tem- 
perature, Although the data is not extremely sensitive to time and 
inlet temperature, provisions will be left in tlie program for time- 
and temperature -dependent data. The above conclusion is based 
on the following: 

a. Pulse shape is not important. For 400-second burns the 
error in not modeling pulse shape is between 0. 005 and 
0. 09 percent. 

b. Thrust varies very little over time. At one hour, a 0. 04 per- 
cent error is introduced by not modeling data as time- 
dependent, At twelve hours, the error only grows to 

0.35 percent. 

c. Inlet temperature variations are within data spread. No 
trends could be seen with the data taken at different tem- 
perature during the performance mapping tests. 
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Table A-1. Repeatability 


To test tlie repeatability of the thrusters used for orbital coutrol, data were 
taken at the beginning and the end of the performance mapping test and were 
compared. At these times the pressure and temperature were the same 
(364 PSIA, 70°F). 


THRUSTER 

NUMBER 

5 

6 

7 

8 

13 

14 

18 

16 

THRUST AT 
BEGINNING 
OF TEST 

0.1225 

0.1240 

0.1225 

0.1245 

0.1200 

0.1265 

0.1260 

■ 

thrust AT 
END OF TEST 

■ 

0.1260 

0.1295 

0.1255 

0,1190 

0.1275 

0.1260 

0.1215 

% change 

1.7% 

1.6% 

4.4% 

0.8% 

0.8% 

0.8% 

0% 

0.4% 


NOTE: HIGHEST PERCENTAGE CHANGE WAS 4.4 PERCENT, 
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Table A-2, Individual Tbruster Modeling 


The differences between the individual thruster test and the perfox'mance map- 
ping test were used to develop the error that could be expected if each orbital 
thruster is modeled. Both high and low pressure points were taken (360 PSIA 
and 160 PSIA); all data were at 70®F. 


THRUSTER 

NUMBER 

6 

6 

7 



8 

13 

14 

15 

16 

INDIVIDUAL 

THRUSTER 

TEST 

160 PSIA 

0,0595 

0.0B30 

0.0615 

0.0615 

0.0590 

0.0630 

0.0615 

0.0595 

360 PSIA 

0.1205 

0.1240 

0.1230 

0.1235 

0.1190 

0.1265 

0.1240 

0.1200 

PERFORMANCE 
MAPPING TEST 

160 PSIA 

0.0625 

0.06/.5 

0.0625 

O.OB25 

0.0595 

0.0635 

0.0615 

0.0600 

360 PSIA 

0.1235 

0.1250 

0.1255 

0.1250 

0.1195 

0.1265 

0.1245 

0.1200 

% CHANGE 

160 PSIA 

4.6% 

0.8% 

1.6% 

1.6% 

0,8% 

0.8% 

0% 

Q.8% 

360 PSIA 

2.5% 

0.8% 

2,0% 

1.2% 

0,ii% 

0% 

0,4% 

0% 

AVERAGE 

3,55% 

0.8% 

0.9% 

1 

1.4% 

0.6% 

0,4% 

0.2% 

0.4% 


NOTE: HIGHEST PERCENTAGE CHANGE WAS 4.6 PERCENT, THIS IS VERY CLOSE TO THE HIGHEST PERCENTAGE 

CHANGE IN TABLE A-1. 
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Table A-3. Average Thruster Modeling 


The' average values of thrust were compared to the individual thruster test data 
and the performance mapping test data to develop an error that could be expected 
if an average thruster value is used. Both hi^ and low pressure points were 
taken (360 PSIA and 160 PSIA.); all data were at 70°F, 


THRUSTER 

NUMBER 

5 

6 

7 

8 

13 

14 

15 

16 

AVERAGE 

VALUE 

160 PSIA 

0.0615 

0.0615 

0,0615 

0.0B15 

0.0615 

0.0515 

0.0615 

0.0615 

360 PSIA 

0.1215 

0.1215 

0.121S 

0.1215 

0.1215 

0.1215 

0.121B 

0.1215 

INDIVIDUAL 
THRUSTER 
TEST % 
CHANGE 

160 PSIA 

3,3% 

2.4% 

0% 

0% 

4.1% 

2.4% 

0% 

3.3% 

3BPPSIA 

0.8% 

2.1% 

1.3% 

1.7% 

2.1% 

4,0% 

2.1% 

1 .3% 

PERFORMANCE 
MAPPING TEST 
% CHANCE 

160 PSIA 

1.6% 

1.6% 

1.6% 

1.6% 

3.3% 

3.2% 

0% 

1.3% 

360 PSIA 

2.6% 

2.8% 

3.4% 

2.8% 

1.7% 

4.0% 

2.5% 

1.3% 

% CHANGE 

160 PSIA 

2.45% 

2.0% 

0.8% 

0.8% 

3.7% 

2.8% 

0% 

2.3% 

360 PSIA 

1.7% 

2.4B% 

2.35% 

2.25% 

1 .9% 

4.0% 

2.3% 

1.3% 

AVERAGE 

2.075% 

2.225% 

1.575% 

1.525% 

2.8% 

3.4% 

1 .1 5% 

1.8% 


NOTE: HIGHEST PERCENTAGE CHANGE WAS 4.1 PERCENT. THIS IS CLOSE TO THE HIGHEST PERCENTAGE 

CHANGE IN TABLE A-1. 
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Table A-4. Comparison of Individual Thruster Modeling 
to Average Thruster Modeling 


THRUSTER 

NUMBER 

5 

6 

7 

8 

13 

’4 

15 

16 

ERROR IN 
INDIVIDUAL 
THRUSTER 
MODELING 

3,55% 

0.8% 

0.9% 

1.4% 

0,6% 

0.<1% 

0.2% 

0.4% 

ERROR IN 
AVERAGE 
THRUSTER 
MODELING 

2.075% 

2.225% 

1.575% 

1.525% 

2.8% 

3,4% 

1.15% 

1.8% 

IMPROVEMENT 
IF INDIVIDUAL 
THRUSTER 
MODELING IS 
USED 

-1,475% 

1.425% 

0.675% 

0.125% 

2.2% 

3.0% 

0.95% 

1.4% 


NOTE: 


THE AVERAGE IMPROVEWIENT WAS 1 .0375 PERCENT WHEN USING INDIVIDUAL THRUSTERS RATHER 
THAN AN AVERAGE THRUSTER. 






































Table A-5. Error in Burn Time Due to Worst Cases 


Taking values from the worst individual thrusters and putting them into the 
rocket equation at both high and low pressure, the error in burn time due to 
the error in ISP and thrust is computed. In the roclcet equation the space^ 
craft's weight was 3000 poxmds with a AV of 0.533 feet per second (typical 
AV maneuver to be expected in stationlceeping maneuvers). 

t ^ WT ISP / ^ -AV \ 

B T \ ® g ISP / 

where t = the burn time in seconds 
B 

WT = the weight of the spacecraft in pounds 
ISP = the ISP either worst case or average in seconds 
T = the thrust either worst case or average in pounds 
AV = the change in velocity required in feet per second 



1B0 PS1A 

360 PSIA 


T 

ISP 



ISP 

D 

AVERAGE VALUE 

0.0615 

217 

808.106 

0.1215 

223.5 

409.026 

WORST CASE 

0.0595 

208 

835.269 

0.1175 

214.5 

422.95 

DIFFERENCE BETWEEN 
AVERAGE AND WORST 
CASE 



27.163 



13.924 

ERROR IN BURN TIME 



3.361% 



3.404% 


NOTE; THIS TABLE SHOWS THAT THE LARGEST ERROR IN COMPUTED BURN TIME USING AN AVERAGE 
THRUSTER WOULD BE APPROXIMATELY 3.5 PERCENT. 
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Table A-6. Final Calibration Values 


The calibration values listed below describe the deviations of eae^ individual 
thruster from the model. The effect of attitude tliruster firings during an 
orbital maneuver is not included. For initial values of the calibration factor, 
this effect is ne^ected. 

The calibration factor talces the following form; 


X 


(1 + C^) 


ave 


where T = the thrust for thruster x 

X 

C* = the calibration f actor /l 00 
f 

"^ave ~ average thrust from the homogeneous model 

The initial value is developed from the test data. After orbital maneuvers 
new calibration factors based on actual in flight thruster performance should 
be developed. The effect of attitude thrusting (approximately 1. 5 percent) 
should also be added to each calibration factor. 


THRUSTER 

NUMBER 

CALISRATiON 

FACTOR 

5 

ons4 

6 

0.78 

7 

1.4 

8 

0.084 

13 

-0.77 

14 

0,084 

15 

-1.2 

16 

-0.34 
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THRUST (-.BSl 


The following homogeneous model of the hydrazine thruster onboard ATS-6 
was developed from the test data; 



150 


35U 

PRESSURE IPSII 


150 


350 

PRESSURE (PSIl 


Figure A-l(a). As a Function of Pressure 




TIME (SEC) TIME (SEC) 


Figure A-l(b). As a Function of Time 
Figure A-1. Final Thruster Model at 70^F 
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